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We study the growth of an JV-component (including N = 1) order parameter when a system 
with a Lifshitz point is quenched from the homogeneous disordered state to the ordered states. 
We study the scaling behaviours of the structure factors for both the non-conserved and conserved 
order parameters in the long time limit after a quench through the Lifshitz point by using the 
large- N and renormalisation group (RG) methods. We construct the analogues of Allen-Cahn and 
Lifshitz-Slyzov growth laws for nonconserved and conserved order parameters which agree with 
our RG results for N — 1. By extending our large- iV methods to the anisotropic Lifshitz point 
we show that the anisotropy is relevant for the growth of the nonconserved order parameter, but 
irrelevant for the conserved case in the large- N limit. We discuss the effects of mode coupling terms 
on scaling of the structure factor for the conserved order parameter case. We also consider the 
ordering dynamics after a quench through an off-Lifshitz point. In a large- N set up we calculate the 
form of the structure factors of a non-conserved order parameter when a system is quenched from 
the homogeneous paramegnetic to the modulated phase. We show that after a quench from the 
paramegnetic to the modulated phase the form of the structure factor violates the standard form of 
dynamical scaling. Our results compare favourably with the available numerical results. 
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I. INTRODUCTION 

The problem of 'phase-ordering kinetics' is the growth 
of order following a temperature quench from the high 
temperature disordered state to the low temperature 
state (for a general review see [1] ) . These systems serve as 
examples of typical non-equilibrium systems. The struc- 
ture factor exhibits temporal scaling after a quench from 
the high temperature phase to the low temperature phase 
of an Ising order parameter. The dynamics depends upon 
whether the order parameter is conserved or not. Fol- 
lowing the success of the large- TV type calculation in the 
critical phenomena (see, for example [2]), the large- TV 
method has been extended to the problem of quench also 
[1]. In the long time limit, this method yields same dy- 
namic exponent as the renormalisation group (RG) ar- 
guments, when the dynamics is non-conserved, but when 
the dynamics is conserved, it yields two time scales - one 
corresponding to the dynamic exponent as given by the 
RG arguments, and a second one, differing from the first 
one by a factor of log(t) where t is the time elapsed af- 
ter quench. However, subsequently, numerical results for 
finite N as well as analytical calculations in a 1 /N expan- 
sion for large- N show that this temporal multiscaling is 
an essentially an artefact of the N — > oo limit [1,3]. For 
finite N one recovers simple scaling. 

In this paper we consider the ordering dynamics of 
non-conserved (NCOP) and conserved (COP) order pa- 
rameters after a temperature quench in a system with 
a Lifshitz point by using large-iV and RG methods. 
We obtain several interesting results: (i) The relevance 



of the anisotropy for the ordering dynamics of NCOP 
but its irrelevance in COP for large-A, (ii) RG analy- 
ses for arbitrary N for both NCOP and COP for the 
isotropic Lifshitz point. Our N = 1 RG results agree 
with our analogues of Allen-Cahn equation for NCOP 
and Lifshitz-Slyzov result for COP, (iii) The relevance 
of mode-coupling terms (torque) for COP with N = 3 
and (iv) the form of the structure factor after a quench 
through an off-Lifshitz point (line EF; see Fig.(l)) in a 
large- N approximation from the paramagnetic phase to 
the modulated phase: We find that the standard form for 
the structure factor exhibiting dynamical scaling breaks 
down. The organisation of the paper is as follows: In 
Sec. II we set up the equations of motion and discuss the 
general scaling form for the structure factor and the two- 
time correlator. In Sec.IIIA we calculate the dynamic 
exponent z and the two-time exponent A after a quench 
through the Lifshitz point for the NCOP by using large- 
TV as well as RG methods. We discuss our analogue of 
the Allen-Cahn result for N — 1. We also elucidate the 
relevance of anisotropy on the ordering dynamics in a 
large- TV set up. In Sec.IIIB we calculate z for the COP 
after a quench through the Lifshitz point by using large- 
TV and RG methods for the isotropic case. We set up 
our analogue of the Lifshitz-Slyzov law heuristically and 
show that it matches with our RG result for N = 1. We 
then show that the anisotropy is irrelevant for the COP 
in the large- iV limit. We also discuss the relevance of 
mode-coupling terms in the equations of motion of COP 
in a simple set up. In Sec. IV, we use the large- N method 
to calculate the structure factor after a quench from the 
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paramagnetic to the modulated phase. In Sec.V we con- 
clude. 



II. EQUATION OF MOTION FOR OUR SYSTEM 

The Landau-Ginzburg free energy for a system with 
(d — m) Lifshitz point is [4,5] 

F = hj ^N> 2 + c ii( v ii^) 2 + c i( v ^) 2 

+ £(V 2 </>) 2 ] +u J d d x(cf > 2 ) 2 ,4> 2 = cj>„cf>„; v=l, N, (1) 

where x = (xj_,X||) is divided into m perpendicular and 
(d — m) parallel components, (d — m) Lifshitz point is 
known as the m-axial Lifshitz point; the special case for 
m = d, i.e., (d — d) Lifshitz point is known as isotropic 
Lifshitz point. The field <j> can be a scalar or an TV- 
component vector (N = 1 corresponds to the scalar field 
case). When both C|| and c± are positive, the ordered 
phase is spatially uniform, but when c±_ < 0, the sys- 
tem can lower energy by creating spatially modulated 
structures with wavevectors | k \ — \ c± \ /2D. The 



point r = 0, cj_ = is a Lifshitz point. There are three 
phases in the phase diagram: The paramagnetic disor- 
dered phase, the low temperature ferromagnetic phase 
and the modulated phase. A schematic phase diagram 
for a vector order parameter of a system with a Lif- 
shitz point is shown in Fig.(l). An example of Lifshitz 
point is the system of mixtures of homopolymers and di- 
block polymers [6]. It also arises in metamagnets with 
appropriate choice of exchanges. A well-known model 
is the Anisotropic Next Nearest Neighbour-interaction 
Ising (ANNNI) model. There have been a great deal of 
work on the static critical proprties at the Lifshitz point; 
recent references include Frisch, Kimball and Binder [7] , 
Diehl, Shpot and Zia [8] , Shpot and Dichl [9], Leite [10], 
Plcimling and Henkel [11]. There have been few studies 
on the dynamical properties at the Lifshitz point; see, 
e.g., Hubcr [12], Folk and Selke [13], Selke [14], and 
Selke [15] . Growth of order after quenches in the ANNNI 
model has been studied by Kaski, Ala-Nissila and Gun- 
ton [16]. Our results are complementary to these works. 
Here, we study quench along the lines CD (cj_ = and 
r changes sign from positive to negative) and EF (cj_ 
changes sign from positive to negative) . For the (d — d) 
Lifshitz point, all C| are zero. 




FIG. 1. A schematic phase diagram for the Landau free energy Eq.(l) for a vector order parameter. Here the high-temprature 
paramagnetic phase (P), the low-temperature ferromagnetic phase (FM) and the low-temperature modulated phase (M) meet at 
the Lifshitz point (LP). Arrows AB, CD and EF refer to quench-paths. The phase boundaries and quench paths are schematic 
only. 



The equation of motion for NCOP is §£ = -T^ 
where as, for COP, L above is to be replaced by TV 2 . 
Since there is no time translation invariance of the 
system, the two-time correlation function C(k, t, t') = 
(</>k(£)0-k(i')) wm be a function of t and t' separately, 
unlike e.g. in dynamic critical phenomena where it is a 



function of t — t'. In the asymptotically long-time limit 
C(k, t, t') typically exhibits a scaling form 

C(M,f) = L(t) d (L(t')/L(t)) x g(kL(t)) 

= t d ' z (t'/t) x ' z g{kL{t)), (2) 

where L(t) ~ t x / z is the characteristic length scale of the 
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system at time t, z is the dynamic exponent, A is the 
two-time exponent and /(0) is a constant. This scaling 
form is expected to hold when both t, t' » t Q where 
t a is some microscopic timescale and L(t) ^> L(t'). The 
structure factor or the equaltime correlation function is 
given by S(k,t) = (0(k, t)<j>{— k, t)). In the long time 
limit, S(k,t) is expected to exhibit a scaling form 



S{k,t) ~ t d l z g{kt x l z ), 



(3) 



where z is the dynamic exponent. We calculate z and 
A in the long time limit given the equations of motion 
above. In the large-iV limit, the scalar field (j> should be 
replaced by <p v with v = 1,..,N. Expressions (2) and (3) 
have been shown to hold good after quenches where the 
initial (high temperature) phase is random and the final 
(low temperature) phase is uniform. However, in a sys- 
tem with a Lifshitz point the following possibilities arise: 
a)A quench from the paramagnetic to the ferromagnetic 
phase through an off-Lifshitz point (along the line AB in 
Fig.(l): The ordering dynamics here is same as that of an 
Ising order parameter which has been discussed in details 
in the literature, b)A quench from the paramagnetic to 
the ferromagnetic through the Lifshitz point (along the 
line CD), and c)A quench from the paramagnetic phase 
to the modulated phase (along the line EF). Expressions 
(2) and (3) are expected to hold for case b); however it 
is not apriori clear if they hold true for case c) also. We 
investigate that here. 



III. ORDERING DYNAMICS AFTER A QUENCH 
THROUGH THE LIFSHITZ POINT 

A. Nonconserved order parameter 

We first calculate the structure factor in the long time 
limit after a quench through the Lifshitz point (along the 
line CD in Fig.(l) in the large- TV approximation. In this 
limit, the equation of motion for the (d—d) Lifshitz point 
becomes (after scaling away r) 



dt 



at T74/ , 

- = -V (pa + ' 



~ -jy(<t>v<t>i<)<i>a 



(4) 



The initial conditions are given as (<fo(k, 0)) = 
0, (<fo(k, 0)(f)j(— k, 0)) = ASij, where k is a wavevector. 
We solve Eq.(4) in a self-consistent one-loop approxima- 
tion which is equivalent to approximate Eq.(4) by (exact 
in the N — > oo limit) 



f^-V 4 + a(t)0, 

This amounts to replace 4> u <j) v /N by its average < 
in the equation of motion. Here </> stands for any of the N 
components of <p a . This equation can also be derived in 



(5) 



> 



a diagrammatic method [17]. Here a(t) has to be solved 
selfconsistently. The formal solution of Eq.(5) is given by 



hip) = fe (O)exp[-fc 4 t + 6(i)], 



(6) 



where b(t) = f*dt'a(t') implying a(t) = = 1 - 
A^ fc cxp[— 2k 4 t + 2b(t)]. The prescribed initial condi- 
tion has been used to eliminate ^(0). We claim that 
a(t) is negligible in the large t limit (which we will jus- 
tify a posteriori). Hence, on putting a — db/dt = we 
get 



1 = A J d d kexp[2b - 2tk 4 } = Cexp(2&)r 



d/4 



(7) 



where C is a constant. We scale out the t dependence by 
using a scaled variable y = tk 4 and absorbed numerical 
factors including A within t. Hence, b = (d/8) lni, which 
immediately gives 



S(k,t) =t d/4 exp(-2/c 4 t), 



(8) 



which is of the conjectured scaling form. It is easy to see 
that the dynamic exponent z=4. Now with the obtained 
value of b, a = ^| <~ \ — > as t — > oo, which is the 
a posteriori justification of the assumption a(t) — > as 
t — > oo. Our analysis immediately tells us that in the 
long time limit, 

Mt) = M0)t d/8 exp[-fc 4 i]. (9) 
This gives, for the unequal time correlation function 

C(fc,t,i') = (0kW^-k(i')>=i' i/4 (y) 1 cxp[-k 4 (t + t% 

(10) 

thus giving X — d/2 same as for an Ising order parameter 
in the N — > oo limit [1]. It will be very interesting to 
see whether A for this system remains the same as that 
an Ising order parameter even after inclusions of 0(1/N) 
corrections. 

For NCOP, however z cannot be calculated exactly 
(unlike COP; see below) through RG arguments, but 
these arguments lead to the existence of another length 
scale in the problem [1,18]. In the Langevin equation for 
NCOP we add a white noise representing thermal noise. 



dt 



o<p a 



(11) 



with (Ca(ki,ti)&(ka,*2)> - 2rT<5 Q ,„<5(ii-i 2 )<J(ki + k 2 ) 
such that the canonical distribution be recovered in the 
equilibrium. Let under rescaling k = k'/b,t' = b z t 



a (k,t) = <f> a (k'/b,b z t')=bX0' a (k',t'). 



(12) 



Hence the structure factor scales as S(k, t) = 
b 2 x -d t d/z g ( kt i/zy The choice of x = d j2 ensures that 
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< cj>'(k',t')<f)'(-k',t') >= S'(k',t') = S(k,t). At the T = 
fixed point, the Hamiltonian is not scale invariant, but 
has a non-zero scaling dimension y: F({b x (j)}) = b v F{(j)). 
After rescaling the Langevin equation reads 

^--^^#^ = -^ = C(^0/r' d3) 

where the new noise C a (^'^') = b x ~ y ( a (k,t) 
and the noise correlator (( a {k', t')'C, v {— fc', 0)) = 
b d - z - 2y 2TT5(t'). Thus we find, using X = d/2, 



i/r' = b d -y- z h(b)(i/T), 

T' = b~ v T. 



(14) 
(15) 



where h (b) is the contribution due to the elimination of 
the small scales. The scale dependence of temperature T 
appears essentially as the coarse-grained Langevin equa- 
tion has been rewritten in the form (13). The Hamilto- 
nian is dominated by the (V 2 </>) 2 term. One can construct 
the Porod law [1] for this problem which tells us 



S(M) 



l 



L N k d+N- 

The energy density e is given by 
e~ (I V 2 </> | 2 ) 

fc 4 (0 k </)_ k ) 

d d k- 



(16) 



' L N k d+N 

=>e~ L~ 4 for N > 4; 
e ~ L -N ( i-N for N < 4 



(17) 



where £ = 1/ fc ma a; is the UV cutoff. It is obvious that for 
large iV, y = d—4. For TV < 4, the dominant contribution 
comes from L (which is the lower cut off or the system 
size) Hence, for N = 1, y = d — 1; N = 2, y = d — 2 
and for TV = 3, y = d — 3. For N = 4, y = d — 4 and in 
addition to that there is an additional logarithm factor. 

Thus for N — 1 we substitute y = d — 1 and obtain 
-=7 = fe 1_z ft,(6)i. Under coarse graining in the free the- 
ory r change only under rescaling: since z = 4 in the 
free theory V ~ b 3 T. The nonlinear term becomes im- 
portant when L(t) 3> w the interface thickness; in that 
regime T reaches its fixed point value T* which can be 
determined by matching with the value in the free the- 
ory (when L(t) < w) [1]. This gives Y* = w 3 T (oj is 
measured in units of the lattice spacing). We now argue 
this (which in turn suggests z = 4) is consistent with 
the more conventional dynamical scaling at the critical 
point T = T c . At T < T c the domain length L(t) ~ 
(Sw 3 rt/M 2 ) 1 / Z where £ is the surface tension and M 
is the magnetisation (see [1] for the corresponding Ising 
case). Thus for z = 4 we have L(t) - (Sw^i/M 2 ) 1 / 4 . 
Here E, T etc. are functions of T. In the critical region, 



S - <r (d ~ 1} , M 2 - £-(<*-4+'7) 5 r - (i-v-zc where ^ is 
the correlation length, rj is the anomalous dimension and 
z c is the critical dynamical exponent of NCOP. Also the 
interfacial thickness ui <~ C- Substituting we find 



L(t)^C(VC Zc ) 1/4 , 



(18) 



which is consistent with the conventional dynamical scal- 
ing [1]. This provides justification of z = 4 and the ex- 
istence of another length scale lu in the problem. This 
can also be seen in an Allcn-Cahn type approach for 
this model (similar to Ref. [1] for an NCOP in the Ising 

.4 i 

model): We get for a flat equilibrium profile = V'((f)) 
where g is the coordinate normal to the interface. Also, 
noting ^ = — §^§f and by using the equation for the 

equilibrium profile, we find ^ = V 2 V • g where we have 
neglected terms containing higher order derivatives of 4> 
with respect to g for gently curving walls. Here g is a 
unit vector normal to the wall (in the direction for in- 
creasing (f). Since ^ is just the wall velocity v we get 
v = — V 2 V • g which is the analogue of the Allen-Cahn 
equation for this model. V 2 V • g is the Laplacian of the 
mean curvature. Since for a single characteristic length 
L, velocity v - ^ and V 2 V • g - giving L(t) - t 1/4 
for non-conserved scalar order parameter. So we find 
z = 4 same as predicted by the RG analysis. 

For TV = 2,3, i.e., for vector order parameters, the 
Hamiltonian renormalises as a free theory but eventually 
flows to the strong-coupling fixed point characterised by 
(f> 2 = const. , which eliminates one degree of freedom from 
the problem but the remaining ones still rcnormalise as 
free fields suggesting that T changes only due to rescal- 
ing. The analogue of Eq. (18) for a vector order parameter 
is L(t) ~ (psTt/M 2 ) 1 /* . Here p s is the large scale spin 
stiffness and p <~ £~( d ~ 4 ) in the critical region. Thus we 
find, for a vector order parameter 



L(t) ~ C(*/C Zc ) 1/4 , 



(19) 



which is again consistent with the conventional dynami- 
cal scaling. So we see z = 4 for both N = 2 and iV = 3. 
The N = 4 case is interesting due to an extra logarithm: 
L(t) - *V 4 [l + 0(1/ In t)], similar to the N = 2 case of 
non-conserved Landau- Ginzburg theory. However, RG 
cannot detect that. 

It is easy to extend the large- N method to the (d — m) 
Lifshitz point: The equation of motion is 



d(j){k,t) _ 2 



Vi^ + o(t)0, 



(20) 



a(t) being the same as in Eq.(5). Neglecting db/dt (like 
the previous case), we get 

l = exp(26) j d d - m fc||d m fc ± exp(-2fc 2 ^)exp(-2fc 4 _t), 

(21) 
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giving b — [(d — m)/A + m/8] Int. Hence the structure 
factor scales as 

S(k, t) = t (<i-™)/2^/4 exp (_2jfe2 t _ 2k\t). (22) 

We see that there are two dynamic exponents zu = 1/2 
and z± = 1/4 for the parallel and the perpendicular di- 
rections respectively, which follows easily from the form 
of the bare propagator for the equation of motion for 
the (d — m) Lifshitz point is G(k,t) = cxp(— k 2 t — kj_t). 

This can be written in a scaling form g(k^t 1 ^ 2 , k^t 1 / 4 ). 
The existence of two dynamic exponents is clear from 
this form. The nonlinear terms in the long time limit 
only modify the scaling form keeping the dynamic expo- 
nents same (within this self-consistent calculation). The 
results above immediately yield surprising results for A: 
As before we find 

M*) = MQ)t (d ~ m)/4+m/8 cxp[-fc 2 t - kit]. (23) 
Thus the two-time correlator scales as 

/ ./\ (d-m)/i / ,/\ m/8 

C(M,t')=i (d - m)/2+m/4 (y) (j) 

xexp[-k 2 (t + t')-k 4 ± (t + t% (24) 

Thus we identify two X's: A|| = (d— m)/2 and Aj_ = m/2. 
So in presence of anisotropy, like the dynamic exponents, 
one finds two two-time correlation exponents. Again it is 
important to find out how 0(1 /N) corrections may affect 
the effects of anisotropy on the A-exponents. 

Folk and Selke [13] studied the relaxation near the Lif- 
shitz point. They characterised relaxation rate by two 
exponents: a single dynamic exponent z and a second 
exponent x which characterises spatial anisotropy. Our 
Eq.(22) can also be written in terms a single dynamic 
exponent z and a new exponent x: 

S(k, t) = t( d - m )/ 2+m / 4 cxp[-2fcfjt - 2k 4 erp t] 

= t(^™)/«+™/(«) exp[-2fcfi(l + fc£ erp /fc||)] (25) 

with z = 2 and x = 2. Thus our results agree with the 
scaling form at the Lifshitz point. 

B. Conserved order parameter 

Having considered the simpler case of quench of the 
nonconserved order parameter, we now go over to the 
more interesting case of conserved order parameter. We 
first work out the (d — d) Lifshitz point using the large- A'' 
method: In the large-A" limit the equation of motion in 
Fourier space becomes [1] 

^ = -fcW + a{t)]<p. (26) 
The solution of this is 



<f>(k,t) = (f>(k,0)exp[~k 2 b(t) - k% (27) 

As usual, b(t) = f*a(t')dt'. We solve for b self- 
consistently. Using the solution of the equation of mo- 
tion we have a(t) = = 1 - A J k d d kexp[2k 2 b - 2k 6 t]. 
Neglecting a(t) in comparison with 1 (which we jusitify 
aposteriori) , 

1 = &(b/tf 4 J dy V d - x exp[2A(y 2 - y% (28) 

with k = (b/tf^y and A = b 3 / 2 /Vi. Evaluating the 
above integral by using steepest descent method we ob- 
tain b = Ct 1 / 3 (lnt) 2 / 3 , (C is a constant) in the long time 
limit. So the asymptotic form of the structure factor is 

S(fc, t )« t -/««(fc/fcm), &m = (^)V6 j (29) 

where 6(x) = 3{a 2/3 x 2 - x 6 ). Now a(t) = db/dt = 
-(C/3)t- 2 / 3 (lnt) 2 / 3 - [C/^t-^l/iXnt)- 1 ' 3 -» as 
t — > oo. Thus there are two time scales present in S(k, t), 
separated by In t; in other words this form, instead of sim- 
ple scaling, exhibits multiscaling similar to the quench of 
conserved ordinary Ising order parameter [1,19]. Follow- 
ing [1] it is easy to show that the equal time correlation 
function C(x,t) = ((f>(x,t) ■ 4>(0,t)} satisfy the equation 

l^=^C-aW 2 (C+^). (30) 

Equation (30) is similar to that in Ref. [20]. This how- 
ever is not correct to 0(1/A^) as pointed out in Ref. [21]. 
This can be derived by using the Gaussian auxiliary field 
method of Mazenko [22]. For N = oo, the C 3 /N term 
can be ignored and the resulting equation exhibits tem- 
poral multiscaling. For finite N we start with a scal- 
ing form consistent with t 1 ^ 6 growth law, i.e., S(k,t) — 
t d/6g( k t d / e ), a{t) = q 2 n /t 2 ' 3 , and C(r, i) = f(r/t^ 6 ), 
where f(x) is the Fourier transform of g(q) and q m is 
a constant. Substituting these in Eq.(30) gives 

| = -(^ + 12 9 5 -12^), ^) = ^(/ 3 ) q , (31) 

where (/ 3 ) q is the Fourier transform of f(x) . We have 
,g(0) = for a conserve order parameter. We have the 
constraint /(0) = 1, i.e., ^2 q g(q) = 1 to fix q m [1]. For 
large N we find that q m <~ (In A) 1 / 3 which diverges with 
A. This defines a lengthscale L(t) <~ i 1 / 6 . Thus, as in 
the Ising case, our approximate analysis gives consistent 
scaling solutions. So we do not find any temporal multi- 
scaling for finite (but large) N. Fourier transform of g(q) 
gives the scaling function f(x). Recently, in a series of pa- 
pers Castellano et al [21,23,24] discussed preasymptotic 
multiscaling and eventual crossover to simple scaling in 
the asymptotic limit for the COP of the Ising model. 
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They showed that there are competitions between terms 
representing multiscaling and simple scaling respectively; 
and consequently they found a crossover time scale be- 
yond which simple scaling is recovered. Analogous anal- 
ysis on Eq.(30) reveals similar features. 

In a RG set up for COP, T should be replaced by 
Tk 2 in the equation of motion that we used while dis- 
cussing the RG arguments for the nonconserved order 
parameter. After suitable coarse graining and rescaling 
of variables, the equation of motion reads (under rescal- 
ing k' = k/b,t' = b z t) 



b 2 - z+2x - y (l/Tk' 2 ) 



SF 



dt' 



+ C k ,(t>)/(Tk> 2 ) (32) 



while the new noise correlator is given by 

< & (*i)& (*2) >= b 2 -^- 2 y2TTk' 2 S(t[ t' 2 ) (33) 

As usual the temperature scales as T' = b~ y T. Since 
the noise correlator is singular in the k — > limit, no 
new term can be generated by integrating out the short 
wavelength fluctuations which is at least as singular as 
the bare correlator [1]. Hence in the long wavelength 
limit, the bare correlator dominates and consequently we 
claim that 1 /V can change only under rescaling. We write 
the recursion relation for 1/T: 



(I/O = b 



2-z+2 X - 



'(i/r). 



(34) 



Provided 1/T is nonzero at the fixed point, we have 
z = 2 + 2\ — y = d + 2 — y where we used 2\ = d. 
As before we have (from Porod law; see previous Sec.) 
y = d-lioTN = l,y = d-2ioTN = 2,y = d-3 
for N — 3 and y = d — 4 for N > 4. Plugging in these 
values we get z — 3 for scalar order parameter (Lifshitz- 
Slyzov law for this model), z = 4 for TV = 2 component 
case, z = 5 for N — 3 and z = 6 for the large-A^ case. 
For N = 4 there is additional logarithm in the form: 
L(t) ~ i 1/6 [l + 0(1/ In t)}; our RG treatment cannot see 
this extra logarithm. It is surprising that we find z for 
the ordering of non-conserved dynamics of the scalar or- 
der parameter (TV = 1) higher than that in the conserved 
dymamics, as normally conservation law is expected to 
slow down the dynamics. Since RG for NCOP is not ex- 
act it is difficult to interpret this surprising result. Nu- 
merical simulations are required to settle this. 

We now show a simple argument leads to the i 1 / 4 
growth for the N — 1 scalar conserve order parame- 
ter (COP). The diffusion current J <~ — V/x with fj, be- 
ing the chemical potential difference between the two 
species. If there is only one scale L(t) at large times then 
V^i <~ Afi/L where Afi is the change in fi. Near coex- 
istence, the free enery difference is of order (j>A/j, (where 
4> is the magnitude of the order parameter in cither of 
the phases). Since in the nucleation of NCOP, there is a 
critical length scale set by a the surface-energy and the 



free 



;e energy difference A/: L <~ a/A,f ~ a/(f>Afi. Thus 
we get J ~ a<pL 2 . Since J <~ <f>L we have 



dL 
dt 



1 

L 2 ' 



(35) 



This gives L(t) ~ t 1 ^ 3 at long times, which is the ana- 
logue of the Lifshitz-Slyzov law for this model. This pre- 
diction is in accordance with our RG result above. 

We extend a simple argument due to J. Das and 
M. Rao [25] to study the effects of the mode-coupling 
terms on the ordering dynamics. In presence of a mode- 
coupling term, the Eq. of motion with the free energy 
(l)is ' 



d<j) a 
dt 



rv- 



SF 



(36) 



Notice that the precisional term involves V 4 instead of 
V 2 as in the precisional dynamics of Hcisenberg ferro- 
magnet. We write Eq.(36) as (in the spirit of dimen- 
sional analysis) a continuity equation and replace the 
corresponding current by the 'velocity' ^ to obtain 



dL 
~dt 



oM 
' L3 ' 



(37) 



where M 0l a and are the equilibrium magnetisa- 
tion, surface tension and spin mobility respectively. It 
is easy to see that beyond a crossover time given by 
tc{G) ~ (T/M g) 5 <~ 1/G 5 dynamics crosses over from 
z = 5 to z = 4. Thus we expect that in presence of a 
torque dynamic exponent is 4, showing the relevance of 
the torque. 

One can easily extend the large-TV calculation to the 
(d — m) Lifshitz point: The equation of motion is 



dt 



= -(ki + kf,)[k 2 + kl +a(t)]<j>. 



(38) 



In the long wavelength limit the solution of Eq.(38) is 
(retaining terms up to 0(fc 4 )) 



cj)(k,t) = (j)(k,0)exp[k 2 b(t) - kit - k\k\t], 



(39) 



with b(t) = J*a(t')dt'. Neglecting a(t) we obtain 
1 = (b/t) d l 2 J d m x ll d d - m x ± ex.p[2/3(x 2 - x\x\ - a; 4 )]. 
Here k = x(6/t) 1 / 2 ,/3 = b 2 /t. By using steepest decent 
method we obtain 

S(k,t) = t WW'W k ^;k m = (^)V4 5 (40) 

where 9'{k/k m ) another function of k/k m . It is easy to 
see that a(t) indeed goes to as t — > oo justifying our 
assumption. Here also we find multiscaling - there are 
two time scales separated by Int. Presumably, this is 
also an artefact of the N — > oo limit as in the (d — d) 
Lifshitz point case. Interestingly, however we do not find 



G 



separate timescales for || and _L directions unlike the non- 
conserved case. The signature of this is already there 
in the bare propagator: G(k,t) — exp[fc 2 (fc| + k\)t] = 
exp[(fcjjfc^ + kpt] in the limit kj_ -> 0,fc|| -> 0. Thus 
under the scaling k — > k/6, t — > 6 4 t, G does not change 
implying z = 4. Here also nonlincarities determine the 
scaling form keeping z — 4. It is interesting to note that 
for the case of (d — m) Lifshitz point for NCOP we ob- 
tain two distinct z's reflecting the anisotropy. However 
for COP the anisotropy does not show up, i.e., one ob- 
tains only one z = 4 in the lowest order. The reason for 
this can be easily understood from the above analysis. 

IV. SCALING IN A TEMPERATURE QUENCH 
FROM THE PARAMAGNETIC TO THE 
MODULATED PHASE 

In the above we have so far considered temperature 
quench through the Lifshitz point from the paramagnetic 
to the ferromagnetic phase. In these cases the initial 
states are random (no order) and the final ferromagnetic 
states are uniform (i.e., the ordering wavevector is 0). 
However, if the system is quenched from high- to low- 
temperatures through an off-Lifshitz point path there are 
two more cases possible: i)From the paramagnetic to fer- 
romagnetic (along the line AB in Fig.(l)): In this case 
the ordering dynamics in the long wavelength limit is 
same as that of an Ising order parameter which has been 
studied extensively; and ii)from the paramagnetic to the 
modulated phase (along the path EF in Fig.(l)). We ex- 
amine the growth of order in the second case in a large-N 
set up. For simplicity we consider the dynamics of the 
non-conserved order parameter (NCOP) only. 

To calculate the form of the structure factor after 
a quench from the paramagnetic phase to the modu- 
lated phase for the NCOP we start with a free energy 
F = fd d x[-±<j> 2 + i(V|,0) 2 - ^(V^) 2 + f (VW + 
u((j> 2 ) 2 ), (j) 2 — <j) v <j>v, v = 1, N for the modulated phase 
we obtain the equation of motion [we consider the general 
case of a system with (d — m) Lifshitz point] 

d(f> _ 5F 

= \7\(j) - Cl y\(t) - DV\<P + a{t)(j). (41) 

Equation (41), as before, is exact in the N — > oo limit. 
Here, a{t) = 1 — ((j) 2 )- The intial condition is 

(0k(O)0_ k (O)) = A. (42) 
The solution of Eq.(41) is 
«£ k (f) = k (O) ex P [-fcft + Cl k\t - Dk 4 ± t + b(t)], (43) 



with a(t) = By using the initial condition Eq.(42) we 
find for the equal time correlation function or the struc- 
ture factor 

S(k,t) = <0k(t)0_ k (t)> 

= Acxp[-2fc 2 t + 2 Cl k\t - 2Dk\t + 2b(t)}. (44) 

Here b(t) has to be evaluated self-consistently. As be- 
fore neglecting a(t) in the long t limit (which we justify 
a posteriori) we find 

1 = A J d d kcxp[-2kp + 2c ± k 2 ± t - 2Dk\t + 26(f)] 

= Acxp(26(t)) J d d - m fc||rf m fc_L exp(-2fc 2 £) 

x exp(2c ± fcit - 2Dk\t). (45) 

The integration over k\\ can be done trivially. For k± we 
use a saddle point approximation in the long time limit: 
We expand F(k±) = c±kj_ — Dk\ about k±_ = = q 
where q Q is the ordering wavevector in the modulated 
phase. In the long time limit we find 

b{t)=F{ qo )Ht^t 1 l% (46) 
This gives S*(k, t) as 

S(k,t) - t :1 ^ ±1 x 

exp[-2fc 2 t - 2(q 2 Q - c ± k 2 ± )t - 2{Dk\ - q 4 )t}. (47) 

Equation (47) clearly violates the standard form Eq.(3) 
for the structure factor as there is a distinct lenght scale 
q^ 1 = \/2D/c_\_. We can identify a dynamic exponent 
z\\ = 2 for the parallel direction by comparing with 
Eq.(3). However for the perpendicular directions dy- 
namic scaling of the standard form is not found. This 
feature manifests more clearly for the (d — d) Lifstitz 
point case where the structure factor in the long time 
limit assumes the form 

5(k, t) ~ t 1 ' 2 cxp[-2( g 2 - c ± k 2 )t - 2(Dk 4 - q 4 )t], (48) 

which clearly does not obey the standard dynamic scal- 
ing form. It is not apriori clear if this is an artefact of 
the large-N method employed here or whether standard 
dynamic scaling is recovered when corrections to 0(1/N) 
arc included. This requires further investigations. It is 
however clear from Eq.(47) that in the long time limit 
along the parallel directions S(k, t) is non-zero for fc|| = 
and along the perpendicular directions S(k, t) is non-zero 
along k± — q a . Thus fcy = and k± = q a are the order- 
ing wavevectors, as expected. It is also clear from Eq.(46) 
that in the large t limit, a(t) — > justifying our assump- 
tion. 

We obtain, from Eq.(43) the two-time correlation func- 
tion for the (d — m) case as 
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exp[-fcf (t + t') -(q 2 - k 2 )(t + t') - (Dk 4 - qj)(t + t')\. 

Thus comparing with the standard form Eq.(2) we can 
identify a two-time exponent for the parallel directions: 
A|| = (d — m)/2. We also detect the existence of a new 
exponent A = 1/4 coming from the perpendicular direc- 
tions which is not of the standard form. Thus, in the 
./V — > oo limit, after a quench from the paramagnetic to 
the modulated phase the standard form of the correlation 
function does not hold. 

Kaski et al studied the scaling properties of the 
anisotropic structure factor after a quench from the para- 
magnetic to the modulated phase [16] in a two dimen- 
sional ANNNI model (d = 2,m = 1 in our notation) 
employing Glauber spin dynamics (i.e., corresponding to 
the NCOP case studied here). There principal findings 
include (i)a dynamic exponent z = 2 for all directions and 
(ii) anisotropic growth rates. These results agree quali- 
tatively with our large- TV results given by Eq.(47) with 
d = 2, m = 1 for which the structure factor is given by 

S(k,t) = tcxp[-2kp-2{ql - c±k\)t - 2(Dk\ - q 4 )t}. 

(50) 

Thus S(0,t) grows as t yielding (in the notation of [16]) 
n s=a 0.5 and for finite wavevectors this growth is modu- 
lated anisotropically, as discussed in Ref. [16]. 

V. CONCLUSION 

We analyse the form of the correlation and structure 
factors after a quench in a system with a Lifshitz point. 
We first consider a quench through the Lifshitz point 
from the paramagnetic to the ferromagnetic phase. We 
studied both the nonconserved and conserved order pa- 
rameters. In particular we obtained the dynamic expo- 
nent z by using large-TV as well as RG calculations. We 
also calculate the two-time exponent A for the noncon- 
served order parameter. We find multiscaling for the con- 
served case in the N — > oo limit and simple scaling for 
the nonconserved case in our one- loop calculations. We 
however show that this multiscaling is an artefact of the 
N — > oo limit and for finite (but large) N simple scaling is 
recovered. We construct Allen-Cahn and Lifshitz-Slyzov 
type scaling laws for the nonconserved and conserved or- 
der parameters respectively. We have also shown that in 
the anisotropic (d— m) case, for a conserved order parme- 
tcr anisotropy is unimportant (irrelevant in a RG sense) 
where as for a non-conserved order parameter it shows up 
in the form of two dynamic and two two-time exponents. 
It is important to see if anisotropy remains unimportant 
in higher order approximations and/or numerical simu- 
lations. We argue that in presence of a torque in the 



Eq. of motion for a conserved order parameter for the 
(d — d) Lifshitz point z reduces from 5 to 4. We have 
also discussed RG arguments in the context of our one- 
loop results. We also show, in the N — > oo limit that the 
standard form of the structure factor exhibiting dynam- 
ical scaling breaks down after a quench from the param- 
agnetic phase to the modulated phase. It is not clear if 
the breakdown of dynamical scaling will persist even af- 
ter the inclusion of corrections of 0(1 /N). This remains 
an open problem. Our results can be easily checked nu- 
merically by extending the cell dynamics methods [26] 
as appropriate for systems with Lifshitz points. Numer- 
ical measurements of higher-order correlation functions 
will help us to understand the validity of the Mazcnko 
method in this problem. Recent numerical studies [25] 
suggest that the Mazenco method may not be very good 
for COP for Ising systems. Numerical simulations are 
required to settle this for systems with Lifshitz points. 
In short, we have analysed the asymptotic form of the 
correlation functions/structure factors after temperature 
quenches in systems with Lifshitz points. We feel our re- 
sults are sufficiently important and interesting to induce 
further detailed studies of these phenomena in systems 
with Lifshitz points. 
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